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ABSTRACT 


A  new  ice  nucleant  aerosol  is  developed  from  the  principles  of  synergism.  Two 

known  ice  nucleants,  silver  iodide  (Agl)  and  bismuth  triiodide  (Bil^)  are  combined 

n  n 

as  2%  AgI-NH4I  X  mole  %  Bil3  NH4I  in  a  water-acetone  solution.  The  resultant 
aerosol  after  combustion  is  tested  for  effectiveness  in  the  Colorado  State  University 

n 

Isothermal  Cloud  Chamber  (ICC),  and  the  optimum  ratio  of  Agl  and  Bil3  is 
determined  to  be  0.5  mole  %  Bits.  The  effectiveness  at  this  optimum  ratio  is  an 
order  of  magnitude  greater  than  Agl  alone,  although  the  temperature  spectrum, 
rate,  and  mechanism  are  similar.  A  subsequent  x-ray  powder  analysis  identifies 
the  optimum  ratio  aerosol  as  the  iodargyrite  form  of  Agl  having  the  closest  match 
to  ice  ever  reported  for  a  hexagonal  crystal.  Unlike  other  mixed  Agl  containing 
aerosols  the  resultant  aerosol  maintained  its  hexagonal  form  of  Agl  at  all  mole 
ratios  of  Bils  tested.  Bismuth  is  at  too  small  of  &  ratio  to  be  identified  through 
the  method  employed  although  it  appears  to  be  fully  incorporated  into  the  Agl 
crystal  as  an  oxide  or  iodo-oxide.  The  plausible  structure  is  looked  at  in  light  of 


epitaxial  and  active  site  theories  of  heterogeneous  ice  nucleation. 
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Chapter  1:  INTRODUCTION 


1.1  Choice  of  an  Ice  Nucleant  Aerosol 

Weather  modification  programs  are  continually  in  need  of  an  ice  nucleant, 
whose  optimum  effectiveness,  in  producing  ice  crystals,  will  occur  at  temperatures 
warmer  than  -20°C  and  approaching  0°C.  At  these  cloud  temperatures  sufficient 
moisture  is  often  available,  although  natural  ice  nuclei  are  seriously  lacking.  While 
the  major  aim  of  research  is  searching  for  an  aerosol  more  effective  at  warmer 
temperatures,  there  is  also  much  to  gain  by  a  more  cost  effective  aerosol.  Silver 
iodide  and  silver  iodide  containing  aerosols  have  long  been  in  operational  use, 
but  the  search  continues  to  modify  silver  iodide  to  be  more  effective  at  warmer 
temperatures  and  more  cost  effective. 

Precipitation  enhancement  programs  must  also  have  an  ice  nucleant  which 
can  be  targeted.  A  nucleant  with  an  unusually  long  time  lag  would  be  of  little 
value  in  most  seeding  operations.  This  lag,  as  measured  by  the  rate  of  ice  crystal 
formation,  is  important  in  determining  the  location  and  type  of  delivery  system, 
as  well  as  the  development  of  ice  crystals  and  the  actual  precipitation  process. 
The  appropriate  ice  nucleant  with  the  desired  rate  to  match  the  intended  target 
will  greatly  aid  the  success  of  a  cloud  seeding  operation. 

Another  factor  in  selecting  of  an  optimum  ice  nucleant  is  its  behavior  in 
particular  environmental  conditions.  The  behavior  of  an  ice  nucleant  is  manifested 
in  its  mechanism  of  ice  nucleation.  The  mechanism  will  reflect  whether  the  aerosol 
is  suitable  for  a  given  cloud  liquid  water  content.  Vapor  concentration  which  is 

temperature  dependent  will  also  greatly  influence  an  aerosol's  behavior. 

I 
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A  final  aim  in  the  choice  of  an  ice  nucleant  aerosol  is  in  understanding  the 
chemistry  of  the  aerosol  catalyzed  ice  nucleation  reaction.  This  is  not  normally 
a  necessary  part  in  the  selection  of  an  aerosol,  but  is  important  in  the  wider 
application  of  cloud  glaciation  operations.  This  is  not  completely  possible  yet. 
but  inference  can  be  drawn  from  chemical  analyses  and  crystallographic  research. 
Knowing  why  a  particular  aerosol  functions  well  as  an  ice  nucleant  might  be  the 
best  lead  to  subsequent  research  to  find  a  better  ice  nucleant. 

1.2  Purpose 

The  objective  of  this  study  was  to  investigate  the  synergistic  interactions  of  two 
ice  nucleants,  silver  iodide  (Agl)  and  bismuth  triiodide  (Bil3).  Synergism  is  the 
joint  action  of  agents  which  when  taken  together  increase  each  others  effectiveness. 
Both  Agl  and  Bil3  were  determined  as  ice  nucleants  almost  forty  years  ago,  but 
little  work  has  been  done  to  promote  Bil3  as  an  operationally  used  cloud  seeding 
agent.  Agl  on  the  other  hand,  has  been  used  in  weather  modification  programs 
since  its  inception. 

Another  purpose  is  to  determine  the  rate  of  ice  crystal  formation  by  the 
resultant  aerosol.  Rates  of  ice  crystal  formation  and  effectiveness  at  several 
different  temperatures  will  also  be  determined.  Such  information  is  a  necessity  to 
make  the  resultant  Agl-Bil3  aerosol  useful  in  a  cloud  seeding  operation. 

A  third  purpose  is  to  better  understand  the  function  of  the  aerosol  as  an 
ice  nucleant.  From  the  effectiveness  and  rate  of  nucleation,  a  mechanism  of  ice 
nucleation  will  be  determined  with  some  reliance  on  chemical  kinetics. 

Beyond  this  macro  study  of  the  mechanism,  the  microstructure  of  the  Agl- 
Bil3  resultant  aerosol  will  be  examined.  Using  an  x-ray  powder  pattern  analysis, 
the  structure  will  be  determined,  and  a  direct  correlation  will  be  drawn  between 
the  crystal  structure  and  the  ice  nucleating  effectiveness  of  varying  ratios  of  the 
Agl-Bil3  hybrid  aerosol.  The  conclusions  drawn  from  this  structure  effectiveness 
study  will  be  used  to  examine  the  epitaxial  (closest  fit)  and  active  site  theories  of 
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Agl  heterogeneous  nucleation,  shedding  some  light  on  these  theories  and  creating 
new  questions  for  further  research. 


\ 
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Chapter  2:  HISTORICAL  PERSPECTIVE 


2.1  History  and  Nucleation  Theory  of  Agl  Containing  Aerosols 

2.1.1  Early  work  with  Agl  and  Bil3 

Vonnegut  (1947)  has  generally  been  credited  with  the  discovery  of  silver  iodide 
(Agl)  as  an  ice  nucleant,  although  he  has  given  credit  to  M.  R.  Block  (Vonnegut, 
1969).  Albeit,  Vonnegut  was  responsible  for  the  pioneering  work  with  Agl  in 
the  weather  modification  field;  he  developed  the  first  solution  combustion  system 
for  the  generation  of  Agl  aerosols.  The  Agl-ammonium  iodide  (NH«I)-acetone- 
water  solution  was  quick-quenched  to  produce  a  submicron  Agl  aerosol  (Vonnegut, 
1949).  The  NH41  was  used  to  solubilize  Agl  in  order  to  prepare  pure  Agl.  The 
NH4I,  used  as  solubilizing  agent,  is  assumed  to  decompose  during  combustion 
(St.  Amand,  1971)  with  no  detectable  incorporation  into  the  Agl  crystal  (Von¬ 
negut,  1949).  Other  solution  combustion  methods  of  Agl  aerosol  production  were 
subsequently  tried  and  used  operationally.  The  2AgI-NaI  acetone  system  was 
commonly  used  during  the  1960’s,  but  neither  this  nor  any  other  acetone  system 
was  as  effective  as  the  original  method  of  Vonnegut's  (Garvey,  1975).  NH4I  as 
the  solubilizing  agent  for  Agl  was  reintroduced  in  1970  (personal  communication 
with  William  Finnegan). 

Sano  and  Fukuta  (1956)  first  discovered  bismuth  triodide  ( Bi I3 )  as  an  ice 
nucleant.  Bil3  was  later  shown  to  have  a  slightly  warmer  threshold  temperature 
than  Agl  as  an  ice  nucleant  in  freezing  bulk  water.  Only  lead  iodide  (Pbl) 
proved  to  be  a  more  effective  nucleant  than  Bi I3  (Fukuta,  1958).  Except  for 
Fukuta’s  theoretical  considerations  explaining  the  nucleating  ability  of  Bil3  . 
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research  operations  have  not  used  BH3  since  1958.  The  theorized  reasons  behind 
the  ice  nucleating  ability  of  Bil3  are  similar  to  Agl.  This  theory  and  the  most 
generally  accepted  theories  of  Agl  nucleation  will  be  discussed  in  the  next  section 
followed  by  a  section  on  the  characterization  of  a  nucleant. 

2.1.2  Heterogeneous  Ice  Nucleation  Theory 

Vonnegut’s  initial  decision  to  use  Agl  as  an  ice  nucleant  was  based  on  the 
similarity  of  the  Agl  crystal  structure  to  that  of  ice.  Such  a  nucleant  might 
better  promote  the  epitaxial  growth  of  ice  crystals  (Vonnegut,  1947).  Crystal 
structure  similarities  are  also  attributed  to  the  nucleating  ability  of  Bil3.  An 
alternate  theory  of  Agl  nucleation  suggests  an  active  site  mechanism.  Evidence 
indicates  impurities  or  defects  on  the  surface  of  a  nucleant  are  responsible  for 
the  adsorption  of  water  in  structural  clusters  (Fletcher,  1960,  1961).  This  theory 
might  also  be  plausible  for  Bil3  nucleation.  However,  it  is  unlikely  that  either 
process  is  solely  responsible  for  the  nucleation  ability  of  an  aerosol.  Any  attempt 
at  bringing  the  structure  of  a  nucleant  closer  to  a  match  with  ice  will,  as  a  matter  of 
consequence,  change  the  surface  properties  of  that  nucleant.  An  intended  change 
of  presumed  active  sites  will  likewise  alter  the  structure  of  a  nucleant.  Therefore, 
the  actual  process  of  heterogeneous  nucleation  is  probably  a  combination  of  these 
two  theories.  An  examination  of  these  theories  in  relation  to  Agl  and  Bil3  will 
follow  in  this  section. 

Just  as  certain  crystals  can  promote  the  phase  transition  of  like  or  similar 
chemical  bodies,  it  was  reasonable  to  assume  this  principle  could  be  applied 
to  other  reactions  involving  a  phase  transition.  Turnbull  and  Vonnegut  (1952) 
related  the  crystallographic  theory  of  crystal  nucleation  catalysis  to  that  of  ice 
nucleation.  In  nucleation  catalysis  of  alloy  production,  the  nucleant  or  catalyst 
crystal  structure  similarity  had  a  direct  relationship  to  the  crystal  structure  of 
the  alloy  produced  and  the  temperature  at  which  the  phase  transition  occurred. 
The  phase  transition  of  water  to  ice  via  Agl  nucleation  showed  the  same  direct 


t 


6 

relationships.  Accordingly,  the  1.5%  disregistry  between  the  lattice  constant  of 
Agl  and  ice,  accounted  for  the  2.5°C  of  supercooling  required  to  freeze  water. 
Turnbull  concluded  that  a  nucleant  with  a  closer  match  to  ice  would  be  an  even 
more  effective  nucleant.  Vonnegut’s  initial  choice  of  Agl  and  its  subsequent  success 
as  a  nucleant  lends  much  credence  to  this  closest  match  theory.  Further  evidence 
in  support  of  this  theory,  was  demonstrated  with  an  ice  nucleant  produced  by  a 
coprecipitating  Agl  and  AgBr  (Chessin  and  Vonnegut,  1971).  The  most  effective 
ratio  of  AgBr  to  Agl  was  40  mole  %  which  showed  50%  less  supercooling  (as 
compared  to  Agl)  to  freeze  water.  X-ray  analysis  of  varying  ratios  of  AgBr  to  Agl 
also  indicated  that  the  heretofore  mentioned  optimum  ratio  of  AgBr  to  Agl  also 
had  the  closest  match  to  ice.  However,  this  match  was  along  the  “a”  axis  only,  and 
the  resultant  crystal  structure  was  no  longer  similar  to  ice.  The  hexagon  crystal 
structure  of  Agl  had  been  changed  to  a  face  centered  cubic  with  the  addition 
of  AgBr.  Chessin  and  Vonnegut  are  careful  to  point  out  that  the  change  in  the 
lattice  structure  of  Agl  by  AgBr  might  also  cause  alterations  in  other  physical 
and  chemical  properties  of  the  nucleant’s  surface  that  may  affect  nucleation. 

The  other  theory,  referred  to  here  as  the  active  site  theory,  was  bom  in  part 
because  of  doubts  about  the  epitaxial  growth  of  ice  and  the  hydrophobic  nature 
of  Agl  (Zettlemoyer,  1961).  Agl,  made  free  of  hygroscopic  impurities  by  react¬ 
ing  pure  silver  with  pure  iodide,  did  not  adsorb  water  and  showed  insignificant 
nucleating  ability  (Corrin,  et  al..  1964).  Doped  Agl,  made  by  purposely  contam¬ 
inating  Agl  with  NOj  showed  a  positive  correlation  between  water  adsorption 
and  increasing  amounts  of  contaminant.  It  was  concluded  Agl’s  effectiveness  as 
a  nucleant  was  due  to  impurities  of  hygroscopic  sites  (Corrin,  1967).  The  idea 
of  contaminated  Agl  is  a  real  occurrence,  since  the  grade  of  Agl  commonly  used 
in  combustion  generation  processes,  is  contaminated  during  its  manufacture  by 
hygroscopic  residual  by-  products,  e.g.  alkali  nitrates  and  halides  (DeMott,  19S2). 
In  addition,  Davis  (1974)  has  shown  the  Agl  aerosol  generated  via  combustion 
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of  the  2%  AgI-NH4I-acetone-water- system  required  a  trace  of  NH4I  for  more 
effective  ice  nucleation. 

Fukuta  (1958)  used  a  crystal  structure  approach  to  explain  the  apparent 
activity  of  Bils -  He  equated  the  “a”  axis  lattice  constant  of  pseudo-  hexagonal 
Bil3  (7.815  A)  to  the  length  of  one  side  of  the  smallest  hexangle  of  ice  (4.5133  x 
v/3  =  7.815  A).  This  might  be  stretching  the  closest  fit  theory  a  bit,  but  it  might 
explain  the  juxtaposition  of  active  sites. 

Clearly,  neither  of  these  theoric--  independently  resolve  all  of  the  questions  of 
heterogeneous  nucleation,  but  likely,  it  is  a  combination  of  both.  The  synergistic 
effect  referred  to  by  Corrin  (1967)  might  be  applicable  to  more  than  a  nucleant- 
contaminant  interaction.  It  seems  plausible  to  look  at  the  synergistic  effect  of 
combining  two  independent  ice  nucleants.  This  study  may  in  turn  shed  more 
light  on  the  theories  of  heterogeneous  ice  nucleation. 

2.2  Characterization  of  an  Ice  Nucleant 

2.2.1  Methods  of  Determining  an  Aerosol's  Ice  Nucleation  Theory 

Several  methods  of  determining  an  aerosol’s  ice  nucleation  activity  or  effec¬ 
tiveness  have  been  tried.  Unfortunately,  some  methods  only  allow  one  mechanism 
of  nucleation  to  operate,  and  all  methods  fail  to  bridge  the  gap  between  the 
laboratory  and  the  real  atmosphere.  Early  attempts  used  cold  boxes  (Shaeffer. 
1946),  but  the  small  volume  to  surface  area  ratio  accounted  for  severe  diffusional 
loss  of  aerosol  to  the  walls.  In  addition,  temperature  gradients  and  the  distur¬ 
bance  caused  by  the  ice  crystal  counting  technique  make  measurements  totally 
unacceptable. 

Following  cold  boxes,  rapid  expansion  chambers  emerged  on  the  scene  (Warner. 
1957).  WJ1  effects  were  reduced  by  treating  the  walls  with  glycol,  but  the  rapid 
cooling  induced  by  rapid  expansion  created  high  transient  supersaturations.  The 
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rate  (time  dependence)  of  the  nucleation  process  was  also  ignored.  This  gave  little 
meaning  to  the  number  of  ice  crystals  formed  (Fletcher,  1958). 

The  membrane  filter  technique  which  utilizes  membrane  filters  to  capture 
aerosol  particles  (Bigg,  et  al..  1963;  Stevenson,  1968),  is  a  method  that  originally 
only  allowed  a  vapor  ice  nucleation  mechanism  to  function.  The  membrane 
filter  is  placed  in  an  ice  thermal  diffusion  chamber  where  the  supersaturation 
of  water  vapor  is  controlled  by  maintaining  a  specific  temperature  difference 
between  the  filter  surface  and  an  ice  surface  above  the  filter.  Though  widely 
used,  it  also  ignores  the  influence  of  the  membrane  substrate  on  the  ice  nucleation 
process  (Fukuta,  1966;  Schaller  and  Fukuta,  1979).  More  recently,  the  membrane 
filter  technique  has  been  improved  to  allow  for  contact  nucleation  by  placing 
water  droplets  on  a  hydrophobic  membrane  filter.  This  technique  has  been  used 
successfully  in  measuring  active  atmospheric  nuclei  (R.  C.  Schnell,  1979). 

Many  other  techniques  have  also  been  employed  to  test  the  effectiveness  of  an 
aerosol,  but  only  large  cloud  chambers  have  been  routinely  used  to  calibrate  the 
effectiveness  of  aerosol  generators  and  pyrotechnics  used  in  weather  modification 
programs.  Colorado  State  University’s  Cloud  Simulation  and  Aerosol  Laboratory 
(Steele  and  Krebs,  1966;  Garvey,  1975)  has  been  a  leader  in  large  cloud  chamber 
design,  testing,  and  operation.  The  large  volume  minimizes  wall  effects  and  allows 
one  to  deal  wdth  the  large  aerosol  samples  generated  from  actual  field  generators. 
As  discussed  in  section  3.1.1,  the  CSU  Isothermal  Cloud  Chamber  (ICC)  does 
not  suffer  from  the  highly  variable  and  high  transient  supersaturations  of  other 
methods.  The  CSU  ICC  can  be  said  to  represent  a  slow  settling  stratus  cloud  at 
water  saturation. 

2.2.2  Rates  of  Ice  Nucleation 

In  the  past,  the  time  dependence  of  ice  nucleation  was  termed  the  “time  lag". 
This  “time  lag”  was  investigated  by  Warner  and  Newham  (1958).  They  found 
the  rate  of  appearance  and  fallout  of  ice  crystals  decreased  exponentially  with 
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time.  Unfortunately,  this  time  dependence  has  often  been  ignored  despite  the 
importance  of  targeting  a  cloud  seeding  operation.  Fletcher  (1958)  understood 
the  importance  of  rate  in  referring  to  the  uselessness  of  effectiveness  without  a 
specification  of  time.  But,  he  incorrectly  interpreted  time  dependence  in  terms  of 
aerosol  particles  activated  as  a  function  of  time  and  temperature  which  assumes 
depositional  behavior.  On  the  other  hand  Isaac  and  Douglas  (1972)  explain  that 
the  same  results  could  have  been  attained  by  contact  nucleation  of  small  aerosols. 

In  all  chemical  reactions,  it  is  understood  that  energy  barriers  induced  either 
physically  or  chemically,  depend  on  a  time  dependent  process.  The  significance  of 
any  reaction  is  highly  questionable  with  an  unacceptable  (either  too  long  or  too 
short)  rate  of  reaction.  Therefore,  any  characterization  of  an  ice  nucleant  aerosol 
is  incomplete  without  a  thorough  understanding  of  the  rate  of  ice  nucleation. 

2.2.3  Ice  Nucleation  Mechanism 

Generally,  four  mechanisms  of  ice  nucleation  are  agreed  upon.  These  are 
defined  here  as  follows: 

vapor  deposition  -  the  deposition  of  water  vapor  directly  to  the  ice  phase  to 
form  an  ice  embryo 

condensation  freezing  -  the  condensation  of  water  vapor  to  the  liquid  phase 
followed  by  freezing  to  form  an  ice  embryo 

condensation  to  droplet  freezing  -  the  complete  condensation  of  water  vapor 
onto  a  nucleus  to  form  a  water  droplet  which  subsequently  freezes 

contact  freezing  -  the  freezing  of  cloud  droplets  upon  direct  contact  by  an  ice 
nucleus 

The  mechanism  of  immersion  freezing  (penetration  of  a  drop  by  an  aerosol 
particle,  followed  by  freezing)  is  not  discussed  here,  because  it  would  be  difficult 
to  separate  from  contact  freezing.  Also,  it  is  unlikely  such  a  mechanism  occurs 
naturally  due  to  surface  tension  effects. 
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While  the  mechanisms  above  and  methods  to  investigate  these  mechanisms 
have  been  discussed  widely,  only  those  methods  which  delineate  nudeation  mech¬ 
anisms  in  competitive  cloud  situations  will  be  discussed.  The  advantage  is  obvious 
from  the  point-of-view  of  weather  modification,  and  affords  the  opportunity  of 
testing  those  aerosols  commonly  used  in  cloud  seeding  operations.  There  are 
four  methods  known  to  have  been  used  for  investigating  in  situ  ice  nucleation 
mechanisms.  The  first  method  analyzed  ice  crystals  obtained  from  an  orographic 
cap  cloud  (Davis  and  Auer,  1972)  and  from  the  CSU  Isothermal  Cloud  Chamber 
(Weikmann,  et  al..  1970)  for  double  structured  frozen  droplet  centered  crystals. 
These  are  associated  with  contact  nucleation  (Auer,  1970).  Davis  and  Auer 
reported  that  increases  in  ice  crystal  concentrations  during  seeding  operations 
with  Agl  (2%  Agl-  NH4I-water-acetone)  -at  -12°C  could  be  explained  entirely 
by  a  collision  nucleation  mechanism.  Wiekmann  et  al.  reported  that  all  aerosols 
tested  are  primarily  contact  nuclei  at  -17°C. 

A  second  method  for  use  in  large  cloud  chambers  used  fluorescent  tracer 
particles  within  cloud  droplets  to  determine  the  percentage  of  ice  crystals  with 
frozen  droplet  centers  (Katz  and  Pilie,  1974;  Katz  and  Mack,  1980).  Only  2-5T 
of  the  crystals  formed  were  the  result  of  contact  nucleation  at  water  saturation. 
Their  work  demonstrated  circular  center  ice  crystals  were  not  unique  to  contact 
nucleation.  Unfortunately,  low  liquid  water  contents  and  an  impractical  nucleant 
(pure  Agl)  reduce  the  significance  of  their  work. 

A  third  method  employed  by  Garvey  and  Davis  (1975)  recognized  the  depen¬ 
dence  of  contact  nucleation  on  particle  size  and  droplet  concentration.  Using 
the  2AgI-NH4I-water-acetone  system,  they  noted  the  differences  in  ice  crystal 
formation  rates  with  changes  in  these  parameters.  The  changes  did  indicate 
contact  nucleation  and  they  obtained  excellent  agreement  between  experimental 
and  theoretical  coagulation  rate  constants. 

The  final  method  used  by  DeMott  (1982)  employed  chemical  kinetic  law  to 
delineate  mechanisms.  He  showed  that  Agl  and  Agl-AgCl  operated  primarily 
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as  contact  nucleants  and  theoretical  calculations  provided  further  proof.  Th is- 
method  is  also  used  here,  to  help  characterize  the  AgI-BiI3-NH4I-  water-acetone 
combustion  aerosol. 
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Chapter  3:  INSTRUMENTATION  AND  PROCEDURES 


3.1  Instrumentation 

3.1.1  Isothermal  Cloud  Chamber 

The  CSU  Isothermal  Cloud  Chamber  (ICC)  has  been  described  in  its  present 
configuration  by  Garvey  (1975).  The  only  modification  is  a  recalibration  of  the 
cloud  inlet  airflow,  an  important  parameter  in  more  accurately  assessing  the  rate 
of  reaction.  This  will  be  discussed  in  section  3.4. 

A  schematic  of  the  ICC  appears  in  Figure  3.1.  A  cloud  is  produced  by  con¬ 
tinuous  atomization  of  distilled  water  with  a  Monaghan  670  ultrasonic  nebulizer. 
As  the  droplets  rise,  they  are  mixed  with  cold  air  and  eventually  escape  out  of 
the  stand  tube  at  the  top  of  the  cloud  chamber. 

The  cloud  density  is  controlled  by  varying  the  rate  of  droplet  introduction 
into  the  ICC.  A  Cambridge  dew  point  hygrometer  monitors  the  cloud  moisture 
content.  The  difference  between  the  saturation  mixing  ratio  corresponding  to  the 
dew  point  temperature  and  the  cloud  temperature  is  the  liquid  water  content 
(LWC).  The  LWC  may  be  varied  between  0.3  and  3.0  gm~3.  The  homogeniety  of 
the  cloud  has  previously  been  examined  with  an  extended  dew  point  probe.  The 
accuracy  of  LWC  measurements  horizontally  and  in  time  has  been  found  to  be 
±  0.15  gm-3  at  1.5  gm"3  LWC.  At  0  5  gm~3  LWC  accuracy  decreases  to  r  0.20 
gm“3.  LWC  accuracies  were  determined  at  -12:C. 

To  maintain  the  cloud  density,  air  must  continually  be  circulated  into  and 
out  of  the  cloud  chamber.  This  air  flow  into  the  chamber  was  measured  at  ear!-, 
temperature  used  in  this  experiment,  and  must  be  taken  into  consideration  in 
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determining  effectiveness  and  rates  of  nucleation.  Measured  air  flow  values  are 
used  to  compute  the  data  in  Table  3.1  and  are  assumed  constant  over  time  for  a 
given  temperature  and  liquid  water  content. 

Temperature  is  monitored  by  an  array  of  thermocouples  and  can  be  maintained 
to  within  ±  0.02°C  at  an  average  chamber  temperature  of  -  16°C.  This  increases 
to  ±0.5°C  at  -20°C.  The  time  response  is  limited  by  the  chart  recorder  to  35 
seconds. 

3.1.2  Vertical  Dilution  Tunnel 

The  vertical  dilution  tunnel  has  also  been  described  in  detail  by  Garvey  (1975). 
The  tunnel  is  operated  in  a  natural  draft  mode  with  the  CSU  standard  test 
generator  below  the  tunnel.  Samples  are  taken  at  the  first  rooftop  platform  with 
a  4.0  liter  syringe.  Dry  (-30°C  dew  point)  air  is  available  to  dilute  the  sample  as 
needed  before  insertion  into  the  cloud  chamber. 

3.2  Data  Collection  Procedures 


3.2.1  Experimental  Procedure 

Much  of  the  procedure  evolved  as  this  study  progressed.  The  first  phase  was 
to  test  and  optimize  the  ratio  of  Bilj  to  Agl  for  best  effectiveness  in  ice  nucleation. 
The  first  solution  tested  was  a  2.0%  AgI-NH4I-20.0  mole%  Bil3-NH4I  in  a  water- 
acetone  solvent.  Each  solution  to  be  tested  is  burned  in  the  CSU  standard  test 
generator  and  the  resultant  aerosol  is  collected  with  a  4.0  liter  syringe.  After 
collection  and  dilution  with  dry  air,  the  resultant  aerosol  is  injected  into  the  ICC 
at  -10°C.  At  3  to  5  minute  intervals,  a  slide  is  pulled  and  examined  in  a  cold  box 
microscope.  Each  slide  is  viewed  at  5  different  slide  locations  with  the  number  of 
crystals  per  viewing  area  counted.  These  five  crystal  counts  are  averaged  to  get 
the  number  of  crystal  per  viewing  area.  This  process  is  continued  until  the  average 
number  of  crystals  per  viewing  area  is  less  than  one  crystal  on  two  successive  slides. 
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The  entire  process  is  repeated  until  the  cumulative  crystal  count  on  two  or  more 
experiments  are  within  10^  of  the  mean.  Subsequent  solutions  with  decreasing 
mole  ratios  of  Bil3  were  tested  in  the  same  manner  at  -10°C. 

After  determining  the  ratio  of  Agl— Bi I3  whose  resultant  aerosol  has  the  highest 
effectiveness  of  ice  nucleation,  this  aerosol  was  again  tested  over  several  temper¬ 
atures  from  -6°C  to  -20°C.  It  was  also  tested  at  -10°C  for  liquid  water  content 
of  0.5  gm-3  in  addition  to  the  standard  LWC  of  1.5  gm'3.  Changes  in  liquid 
water  content  leading  to  changes  in  a  reaction  kinetic  rate  coefficient  is  the 
primary  method  of  distinguishing  between  droplet  dependent  or  vapor  dependent 
mechanisms. 

3.2.2  Collection  of  Sample  for  X-Ray  Analysis 

After  determining  effectiveness  and  rates  for  each  solution,  results  suggested 
the  need  for  an  x-ray  powder  determination  in  an  effort  to  identify  the  resultant 
aerosols  and  to  intercorrelate  the  results.  Aerosol  samples  were  collected  by 
placing  a  nucleopore  filter  attached  to  a  suction  pump  in  the  dilution  tunnel. 
About  1  hour  was  required  to  collect  1  gram  of  aerosol.  The  collected  aerosols 
were  then  sealed  in  vials  and  sent  to  the  Molecular  Structure  Corporation  in 
College  Station,  Texas  to  identify  and  determine  crystal  lattice  spacings. 

3.3  Procedure  of  Data  Analysis 

3.3.1  Effectiveness 

Effectiveness  is  defined  as  the  total  number  of  ice  crystals  which  will  grow 
and  fall  out  per  gram  of  aerosol.  In  the  previous  section,  the  method  of  counting 
crystals  was  discussed.  The  cumulative  number  of  ice  crystals  counted  during  a 
single  experiment  is  related  to  effectiveness  through  the  following  equation. 

Effectiveness  (#Crystals'g  nucleant)  =  A’,r  x  —  x  —  x  —  (3.1) 

A  t  V ,  Rj 
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Nie  =  cumulative  number  of  ice  crystals 

Ae  =  Cross-sectional  area  of  cloud  chamber  (cm  ) 

Av  =  Microscope  viewing  area  (cm  ) 

D,  =  Syringe  dilution  factor 

V,  =  Syringe  volume  (1) 

Rd  =  Wind  tunnel  dilution  rate  (1/min) 

Rg  =  Aerosol  generation  rate  (g/min) 

Ac,  A„  and  V,  are  constants  given  by  8.35  x  103  cm,  2.01  x  10-2  cm,  4.0  liters, 

respectively. 

D,  accounts  for  the  syringe  dilution  needed  to  allow  nucleation  to  proceed 
without  over  depleting  cloud  liquid  water.  At  colder  temperatures  (e.g.  -  16:C, 
three  dilutions  with  dry  air  are  needed.  At  temperatures  warmer  than  -8°C.  no 
dilutions  are  necessary.  Syringe  dilution  factors  employed  are  8.64,  74.7  and  646.0 
for  one,  two  or  three  dilutions  with  dry  air,  respectively. 

is  determined  from  the  vertical  motion  in  the  wind  tunnel  during  aerosol 
generation.  1.0  x  10s  1  min-1  is  used  for  Rj>.  This  value  can  vary  under  variable 
atmospheric  conditions;  however,  the  resultant  inaccuracies  are  minimized  by 
multiple  testing  at  different  times  on  different  days. 

R,  differs  according  to  solutions  used  and  the  generator  nozzle.  This  value  is 
redetermined  periodically  and  whenever  the  nozzle  is  cleaned  or  replaced. 

The  effectiveness  equation  (3.1)  simply  relates  the  laboratory  situation  to  the 
real  atmosphere.  Some  short  falls  of  individual  parameters  are  noted,  and  in 
most  cases,  repeatability  of  the  data  will  limit  errors.  Other  inherent  problems 
in  modeling  the  cloud  chamber  will  be  discussed  in  section  3.5  under  dilution 
mechanisms. 

3.3.2  Rates 

The  rate  of  reaction  or  ice  crystal  production  is  determined  from  the  data  by 
computing  the  percent  production  of  ice  crystals  versus  time.  In  this  way.  rales 
calculated  from  various  runs  can  be  compared  despite  differences  in  total  number 
of  crystals.  Multiple  tests  at  the  same  temperature  and  for  the  same  aerosol  are 
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averaged  to  obtain  average  and  more  representative  plots.  As  already  mentioned.  i 

variations  in  natural  draft  dilutions  can  affect  both  the  concentration  and  sizes  ! 
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of  aerosols.  While  data  multiplicity  can  limit  the  concentration  differences,  par¬ 
ticle  size  can  have  a  subtle  effect  on  the  nucleation  mechanism.  This  limits  the 
sensitivity  of  the  data,  and  at  present,  arithmetic  averaging  appears  appropriate. 


3.3.3  Kinetic  Approach  to  Rates 

In  an  effort  to  better  portray  and  discern  the  rate  data,  a  chemical  kinetics 
approach  is  used.  DeMott  (1983)  showed  the  power  of  this  technique  in  charac¬ 
terizing  a  mixed  Agl-AgCl  aerosol. 

In  a  reaction 


A  +  B  — t+  IC 


(3.2) 


A  =  Aerosol 
B  =  H20 
IC  =  Ice  crystals 

B  is  held  constant  and  in  great  abundance  so  the  reaction  depends  only  on  A. 


The  kinetic  expression  is: 


_d[A j  ,  * 

dt  ~Koi*A‘ 


(3.3) 


[A]  =  Concentration  of  A 

Ko(l,  =  Kinetic  rate  coefficient 
N  =  Order  of  reaction 

In  actuality,  the  reaction  is  second  order  (the  reaction  being  dependent  on 


both  the  concentration  of  A  and  B).  However,  by  holding  B  in  great  abundance, 


the  reaction  depends  only  on  the  concentration  of  A  (N  =  1),  i.e.  pseudo  first 


order.  Integrating  and  solving  for  A: 


'  d’l. 

=  Koi,t 

]A  „ 


the  concentration  of  A  at  any  time  t 
the  initial  concentration  of  A 


(3.-4) 
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Although  [A  is  not  an  observed  parameter,  it  can  be  put  in  terms  of  the 
product  IC,  (one  A  reacts  with  one  B  to  form  one  IC)  i.c. 


\*o  [ICj 

| IC}t  =  the  concentration  of  ice  crystals  at  any  time  t 
[IC]f  =  the  final  concentration  of  ice  crystals 

Substituting  equation  (3.5)  into  (3.4) 

ln[l  -  1^-1  =  ~Koi,t  (3.(1) 

1 1L  J 

Because  this  equation  was  originally  in  terms  of  [A1,  as  the  reaction  proceeds 
[A]  is  depleted.  Therefore,  the  plot  of  equation  (3.6)  will  have  a  negative  slope. 
Taking  the  absolute  value  of  equation  (3.6)  will  give  a  positive  slope.  This  will 
better  represent  the  production  of  ice  crystals  as  the  reaction  proceeds.  Therefore, 
equation  (3.6)  is  plotted  as  {|  ln\l  —  |}  vs  time  (t)  with  the  previously 

determined  rate  data.  K0js  represents  the  slope;  thus,  changes  in  the  rate  are 
easily  observed  as  a  change  in  slope. 


3.3.4  Mechanisms 

The  kinetics  approach  discussed  in  the  previous  section  assumes  a  single  stage 
pseudo  first  order  reaction.  Such  a  reaction  would  result  in  a  straight  line  behavior 
inferring  a  reaction  or  mechanism  with  only  one  rate  determining  step.  In  the  case 
of  contact  nucleation,  contact  is  a  rate  determining  step.  This  rate  determining 
step  is  second  order,  but  the  overall  reaction  is  pseudo  first  order  by  virtue  of 
experimentally  holding  the  water  droplet  concentration  constant.  Equation  (3.7) 
below  represents  this  process. 

A  -  W'D  ^  \VDS  £  1E^  IC  (3.7) 

Wh  ere  A  represents  an  aerosol,  \VD>-  represents  a  nucleated  water  droplet.  IE  is 
an  ice  embryo  or  frozen  droplet,  and  IC  is  the  ice  crystal.  Kj.  K2.  and  K3  are  the 
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rate  constants  associated  with  contact,  freezing,  and  crystallization  from  the  vapor 
phase  respectively.  K0{„  is  the  o%’eral!  rate  constant  experimentally  determined. 
In  contact  nucleation,  WD*  — ♦  IE  — ♦  IC  is  very  fast  so  K2  and  K3  are  large 
and  do  not  contribute  significantly  to  the  observed  rate  constant.  If  however, 
the  rate  constant  associated  with  freezing  were  to  represent  a  slower  process 
and  approached  the  rate  constant  associated  with  the  contact  step,  straight  line 
behavior  would  no  longer  be  observed.  Thus,  straight  line  behavior  will  effectively 
eliminate  condensation  freezing  and  ‘condensation  to  droplet  freezing’  as  possible 
mechanisms.  The  delineation  between  contact  and  deposition  mechanisms  can  be 
related  to  droplet  concentration  and  will  be  discussed  in  the  analysis  of  the  data. 

3.3.5  Dilution  Correction 

In  adjusting  the  cloud  chamber  to  conform  to  atmospheric  conditions,  the 
problems  with  syringe  dilution  have  been  noted.  Colder  cloud  chamber  temper¬ 
atures  will  show  faster  than  normal  rates.  The  exact  effect  is  unknown.  Other 
factors  which  also  may  affect  the  rate  and  the  effectiveness  are  airflow  dilution 
and  wall  diffusion. 

Wall  diffusion  is  the  diffusional  loss  of  aerosol  to  the  walls  of  the  cloud  chamber. 
It  is  very  dependent  on  particle  size  and  thus  difficult  to  evaluate.  Davis  (1974) 
indicated  wall  diffusion  is  small  compared  to  air  flow  diffusion  and  becomes  much 
more  negligible  for  higher  air  flows  and  faster  rates.  Consequently,  wall  diffusion 
is  not  considered  a  major  factor  and  is  therefore  neglected  in  this  study. 

Air  flow  dilution  is  independent  of  particle  size,  but  very  dependent  on  temper¬ 
ature  and  liquid  water  content.  It  can  be  determined  quite  easily  (Davis,  1974). 
The  kinetic  expression  for  the  loss  of  aerosol  due  to  air  flow  dilution  is: 


20 


Integrating  and  solving 


-,n^D  =  vl 


Equation  (3.8)  is  very  similar  to  (3.6).  ^  is  the  slope  due  to  air  flow  dilution. 
Equation  (3.9)  and  (3.4)  can  be  combined  to  determine  the  corrected  kinetic  rate 
coefficient  (Kcor): 


lnn  .■  .1  \AD  ~  [fcobs 

1  A\o  [A,t  v 


(3.10) 


The  corrected  rate  expression  becomes: 


Where 


\A\ 

-\lnhr\cor  =  Kcort 

[A\0 


Kcor  —  f^obs  ,, 


(3.12) 


Thus,  the  data  plotted  with  equation  (3.6)  is  replotted  using  Kcor  as  the  new 
slope.  Values  for  ^  appear  in  Table  3.1. 

Table  3.1:  Air  Flow  Dilution  Rates 


Temperature 

LWC  (girT3) 

^  (min-1) 

-8 

1.5 

.0473 

-10 

1.5 

.0579 

-12 

1.5 

.0651 

-16 

1.5 

.0693 

-20 

1.5 

.0751 

-10 

0.5 

.0297 

Effectiveness  can  also  be  corrected  for  air  flow  dilution  using  the  same  method. 
Rewriting  equation  (3.4)  and  solving  for  A  ( 

[A]  -  \A\0eip K^.t (  (3.13) 
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Digressing  from  the  corrected  kinetic  rate  expression  (3.11)  to  a  differentia)  form 
and  noting  again  that  one  aerosol  is  assumed  to  produce  one  ice  crystal. 


iA-K  A  =il± 

dt  c"r‘-  '  dt 


(3.11) 


Substituting  (3.13)  into  (3.14): 

d'lC' 


dt 


-  =  Kcor[A]0exp'  -  Kobft] 


(31 1 


Integrating  and  solving  for  )IC) 


K 


[■ IC}t  =  ■~~\A\€xp\  —  Kobft] 

n0b> 

As  the  reaction  comes  to  a  close,  t— *oc.  Letting  t— *oc  in  (3.16) 


(3.16) 


K ; 


i  IC}f  =  Tr 


*  obs 


\0 


Rearranging 


l  A\.  =  ^[IC\ 

**  rnr 


(3.17) 


(318) 


Therefore,  multiplied  by  the  observed  effectiveness  ([IC’fofai/g  aerosol) 

will  give  an  effectiveness  corrected  for  air  flow  dilution.  Correcting  for  air  flow 
dilution  gives  a  net  result  of  increasing  the  total  number  of  ice  crystals  since 
K06,  >  K,;0r;  vice  versa,  airflow  dilution  decreases  an  aerosols  actual  effectiveness. 
Syringe  dilution,  wall  diffusion,  and  air  flow  dilution  also  result  in  increasing  the 
rate  by  accelerating  the  depletion  of  aerosol. 


Chapter  4:  DATA 


4.1  Effectiveness  Data 

Effectiveness  or  number  of  ice  crystals  produced  per  gram  of  a  nucleant  aerosol 
is  presented  in  Figures  4.1  and  4.2.  Figure  4.1,  effectiveness  versus  mole  percent 
ratio  of  Bil3  to  Agl,  shows  the  Agl  Bil3  ratio  which  gives  peak  effectiveness.  This 
ratio  is  0.5  mole  %  Bil3  to  1  mole  Agl  in  solution  at  a  mean  effectiveness  of  4.3 
x  1014,  which  is  subsequently  referred  to  as  the  optimum  ratio.  The  small  open 
circles  are  data  points  from  individual  experiments.  In  order  to  better  present 
the  data,  both  ordinate  and  abscissa  are  log]0  scales.  The  effectiveness  of  the 
aerosols  produced  by  combustion  of  a  2%  AgI-NH4I-acetone-water  solution  and 
a  2 %  BiI3-NH4I-acetone-water  solution  are  represented  as  horizontal  dotted  lines. 
These  dotted  lines  show  the  theoretical  left  (Agl  with  no  Bismuth)  and  right  (Bil3 
with  no  Agl)  asymptotes.  The  upper  solid  line  corrects  this  figure  for  airflow 
dilution.  Figure  4.1  is  not  corrected  for  airflow  dilution  since  its  purpose  was 
to  show  the  optimum  ratio  of  Bil3  to  Agl.  Data  points  are  represented  by  open 
circles.  Figure  4.2,  effectiveness  versus  temperature  is  shown  as  the  dotted  line.  A 
temperature  of  -10°C  and  LWC  of  1.5  gm‘s  are  used  for  the  generation  of  all  this 
data.  -10°C  was  chosen  somewhat  arbitrarily  as  a  place  to  start.  Temperatures 
warmer  than  -10CC  might  have  made  it  impossible  to  detect  the  effectiveness  of 
some  aerosols.  Colder  than  -16CC  would  have  caused  many  of  the  better  mole 
ratio  aerosols  to  approach  the  same  effectiveness  value.  A  LWC  of  1.5  gnf3  is 
chosen  accordingly  in  addition  to  providing  a  competitive  cloud  environment  of 
water  droplets  and  vapor. 


Rates  of  ice  nucleation  are  shown  in  Figures  4. 3—4. 7.  Individual  data  point" 
are  mean  aveiage«  of  all  runs  at  that  temperature.  Figure  4.3  and  Figure  4.4  show 
the  percent  of  reaction  completed  versus  time.  The  former  compares  the  rate  of 
ice  nucleation  with  Agl  and  the  ice  nucleation  rate  of  the  optimum  ratio.  Roth 
indicate  similar  rates.  Figure  4.4  demonstrates  rate  as  a  function  of  temperature  . 
LWC  remains  constant.  Kinetic  plots  (equation  (3.6))  of  the  same  data  are  shown 
in  Figures  4.5  and  4.6.  While  it  is  difficult  to  see  any  significant  difference  in 
the  rates  of  Figures  4.3  and  4.4,  their  respective  kinetic  plots  (Figures  4.5  and 
4.6)  provide  greater  resolving  power  for  the  rate  differences.  Figure  4.7  is  a 
kinetic  plot  indicating  rate  as  a  function  of  LWC.  This  plot  will  serve  as  an  aid  in 
determining  the  mechanism  of  nucleation.  The  term  raw  data,  in  figures  4.5-  4.7 
refers  to  the  data  not  having  been  corrected  for  airflow  dilution.  Figures  4.8-4.10 
corresponding  to  Figures  4. 5-4. 7  are  kinetic  plots  corrected  for  airflow  dilution. 

4.3  X-ray  Powder  Analysis  Data 

Three  samples  of  varying  mole  ratios  of  Bil3  to  Agl  wre  chosen  for  x-ray 
powder  analysis.  Samples  were  collected  and  sent  to  the  x-ray  laboratory  a" 
mentioned  in  section  3.4.  The  three  samples  chosen  were  the  aerosols  produced  via 
combustion  of:  the  optimum  ratio  (0.5  mole  percent  B i I3  to  Agl):  Agl  only:  and  a 
10  mole  percent  ratio  of  Bi I3  to  Agl.  From  figure  4.1  these  three  aerosols  represent 
a  peak  effectiveness,  and  left  and  right  extremes  of  effectiveness,  respectively. 

The  data  received  included  d-values  (crystal  lattice  spacing")  with  maximum 
probable  errors;  corresponding  relative  irradiances,  2-Theta  values,  and  peak 
widths  for  each  sample.  Table  4.1  presents  this  data.  Computer  identification 
of  each  aerosol  sample  indicated,  in  all  case".  Agl  IodaTg\  rite,  i.e  the  hexagonal 
crystal  form  of  Agl  as  opposed  to  the  cubic  form.  Analysis  of  this  data  and 
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reduction  of  lattice  spacing*  (d-values) 
next  chapter. 


vV/v'.-’/vV1 


25 


I 

t 


I 

l 


Effectiveness  (♦  ice  crystols/g  nude  ant) 


CSU  Standard  Generator 

2%AgINH«I  /  0.5  mole  V.Bil 
LWC  =  l.5gnr3 
Nolurol  Droft 


0-2  -4  -6-8  -10  -12  -14  -16  -18  -20  -22  -24 

TEMP  (#C) 


Figure  -1.2:  Effectiveness  of  nucleation  versus  temperature  for  the  aerosol  produced 
via  comhusion  of  the  optimum  ratio  solution  of  Agl  and  Bil.v  1  tie 
dashed  line  represents  the  best  fit  to  the  observed  data  poin’1-. 
and  the  upper  solid  line  represents  the  observed  data  corrected  for 
airflow  dilution. 
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Figure  4.3:  Rates  of  ice  crystal  formation  comparing  the  aerosols  produced  via 
combustion  of  Agl-NH^-O.o  mole  Sc  BH3-NH.4I. 
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Figure  4.7:  Raw  kinetic  plots  varying  the  liquid  water  content  in  the  cloud 
chamber 
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Figure  4.8:  Kinetic  plots  of  Figure  4.5  corrected  for  airflow  dilution 


Figure  4.9:  Kinetic  plots  of  Figure  4.6  corrected  for  airflow  diiuti 
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Figure  4.10:  Kinetic  plot  of  Figure  4.7  corrected  for  airflow  dilution 


Table  4.1:  Data  from  x-ray  powder  analysis,  Wavelengths  used  1.54184,  1.54056, 
and  1.54439. 

2 %  AgI-NH4I-0.5  mole  %  Bil3-NH4I 

D  value  I/IMAX 

Dj:3. 96607  ±  .00350  33.60 

D2:3. 74654  ±  .00311  100.00 

D3:2. 29252  ±  .00112  69.10 

2%  AgI-NH4I 

D  value  I/IMAX 

0^4.00237  ±  .00356  36.10 

D2:3.77846  ±  .00317  100.00 

D3:2.30544  ±  .00114  76.80 


2%  AgI-NH4I-10  mole  %  BiI3-NH4I 

D  value  I/IMAX 

Di:4.01007  ±  .00358  47.00 

D2:3. 78373  ±  .00318  100.00 

D3:2. 30681  ±  .00114  83.00 
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Chapter  5:  ANALYSIS 


5.1  Effectiveness 

One  of  the  most  striking  features  of  ice  nucleation  by  the  AgI-BiI3-NH4I 
sy  tern  is  seen  in  figure  4.1.  The  highest  effectiveness  occurs  at  a  ratio  of  1  mole 
Agl  to  0.5  mole%  Bil3.  This  kind  of  specific  mole  ratio  of  one  compound  to 
another  compound  has  also  been  demonstrated  in  other  ice  nucleation  aerosols. 
DeMott  (1982)  produced  a  mixed  Agl-AgCl  aerosol  via  solution  combustion  of  a 
2%  AgI-NH4I-acetone-  water  with  varying  mole  ratios  of  NH4CIO4.  He  attained 
the  highest  effectiveness  with  a  ratio  of  30  moIe9c  NH4CIO4.  Vonnegut  and  Chessin 
(1971)  with  their  Agl-AgBr  aerosol  produced  the  highest  effectiveness  with  a  ratio 
of  40  mole%  AgBr.  However,  the  relative  proportions  of  reactants  Agl  and  Bil3 
is  different  from  previous  Agl  containing  aerosols  in  that  the  mole  concentration 
of  silver  is  200  times  that  of  bismuth.  The  aerosol  of  this  ratio  also  has  an 
effectiveness  an  order  of  magnitude  greater  than  Agl  alone,  which  demonstrates 
the  dramatic  change  such  a  small  portion  of  another  compound  (Bil3)  can  have 
on  the  resulting  ice  nucleant.  The  reasons  for  this  increased  effectiveness  by  the 
optimum  ratio  of  Agl/Bil3  resultant  aerosol  will  become  clear  upon  examining 
the  morphology  of  the  resultant  aerosol  in  Section  5.3. 

Figure  4.2  shows  the  optimum  ratio  approaching  a  limiting  effectiveness  of 

1.1  x  1015  at  -20:C.  This  is  not  surprising,  however  since  all  nucleants  approach 
some  minimum  temperature  at  which  a  limited  effectiveness  is  attained.  Better 
effectiveness  at  warmer  temperatures  is  a  common  goal  for  weather  modification 
programs,  but  -6CC  is  the  warmest  temperature  used  in  this  study.  It  is  clear. 
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however,  that  the  optimum  ratio  Agl  Bi I3  resultant  aerosol  is  not  a  warm  tem¬ 
perature  candidate. 

5.2  Rates  and  Mechanisms 

Raw  rates  of  ice  crystal  formation  are  shown  in  Figures  4.3  and  4.4.  It  is 
immediately  obvious  from  Figure  4.3,  that  the  standard  Agl  aerosol  is  somewhat 
mimicked  by  the  optimum  ratio  Ag]/Bil3  aerosol.  Thus,  the  rate  and  mechanism 
are  probably  similar,  but  the  efficiency  of  ice  crystal  catalysis  (figure  4.1)  is  much 
greater  for  the  optimum  ratio  aerosol.  The  process  of  ice  crystal  formation  is  a 
very  slow  process  taking  30  minutes  or  more  for  completion,  but  as  much  as  90r~ 
of  the  reaction  is  complete  in  less  than  15  minutes. 

Figure  4.4  clearly  shows  that  the  ice  crystal  formation  rates  are  a  function 
of  temperature,  i.e.  the  rate  increases  with  decreasing  temperature.  Also  each 
curve  has  a  similar  exponential  shape.  These  factors  indicate  the  aerosol  probably 
functions  in  the  same  way  at  all  temperatures  and  also  functions  like  the  standard 
Agl  aerosol.  However,  kinetic  plots  are  needed  before  the  true  character  of  these 
rates  can  be  discerned. 

Figure  4.5  and  4.6  are  raw  kinetic  plots  of  the  previous  rate  data.  Indeed,  the 
data  is  easier  to  discern  in  relation  to  rates.  The  ice  formation  rate  of  Agl  and 
the  optimum  ratio  aerosol  in  figure  4.5  are  no  longer  inseparable.  The  rates  are 
similar,  but  are  not  identical.  The  optimum  ratio  aerosol  has  slightly  faster  rate 
of  ice  crystal  formation.  Figure  4.6  demonstrates  the  temperature  dependence 
of  ice  formation  rates  much  more  clearly  than  the  nonkinetic  plots  of  figure  4.4. 
However,  more  important  in  both  figure  4.5  and  4.6  is  the  straight  line  behavior 
of  the  data.  As  mentioned  in  section  3.3.4,  linear  behavior  is  indicative  of  a  single 
stage  pseudo  first  order  reaction.  All  temperatures  indicate  pseudo  first  order 
with  a  singular  mechanism  of  ice  nucleation.  To  determine  the  actual  mechanism, 
the  relationship  between  rate  and  liquid  water  content  must  be  determined. 
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The  liquid  water  content  and  temperature  of  a  supercooled  cloud  parcel  could 
make  a  big  difference  in  the  nucleation  rate  and  mechanism  in  operation.  Of 
course,  the  ICC  allows  one  to  maintain  a  constant  temperature  and  liquid  water 
content  to  better  discern  the  mechanism  present.  As  mentioned  before,  in  a 
constant  volume  saturated  environment  such  as  the  ICC,  a  decrease  in  liquid 
water  will  lead  to  decrease  in  droplet  concentration.  Vapor  concentration  remain^ 
constant  because  it  is  temperature  dependent.  Therefore,  a  change  in  liquid 
water  would  have  no  effect  on  a  vapor  dependent  process.  On  the  other  hand 
a  droplet  process  such  as  contact  nucleation  would  show  a  change  in  rate.  Figure 
4.7  shows  a  direct  relationship  between  liquid  water  content  and  rate  of  ice  crystal 
formation.  Clearly,  at  -10°C  and  all  temperatures  tested,  the  optimum  ratio 
Agl-Bil3  resultant  aerosol  employs  contact  nucleation  in  the  production  of  ice 
crystals.  Contact  nucleation  is  a  psuedo  first  order  (second  order)  single  stage 
reaction.  DeMott’s  (1982)  kinetic  study  of  Agl-AgCl  indicated  many  of  the 
same  results.  However,  at  -  20CC  his  aerosol  demonstrated  a  vapor  dependent 
mechanism  beyond  20  minutes.  If  this  mechanism  exists  for  the  optimum  ratio 
aerosol  it  is  disguised  in  the  noise  since  the  reaction  is  969c  complete  in  18  minutes. 
Otherwise  the  optimum  ratio  Agl-Bil3  aerosol  serves  as  a  contact  nucleator  in  the 
production  of  ice  crystals  with  a  slightly  faster  rate  of  reaction  as  the  standard 
Agl  aerosol. 


5.3  Nucleant  Morphology 


During  the  early  stages  of  this  study,  it  was  suspected  that  Bi I3  underwent 
some  stage  of  oxidation  during  the  quick  quenching  process  after  combustion 
generation.  The  first  real  clue  came  in  the  erratic  data  generated  when  the 
3%  AgI-NH«I-x  mole  9£  Bi I3  solution  remained  idle  for  a  period  of  time.  In 
addition,  a  flask  of  solution  lost  its  characteristic  ured  Kool-Aid"  color,  turning 
brown  and  murky,  after  only  one  day.  The  aerosol  generated  via  combustion  of  a 
0.259c  BiI3-NH«I  solution  was  identified  from  color  and  melting  point  as  B^O.c. 


confirming  BiU  readily  oxidizes.  As  an  ice  nucleant,  Bi203  demonstrated  minima! 
effectiveness  (Figure  4.1).  There  is,  however,  no  indication  the  bismuth  form 
incorporated  into  the  crystal  structure  of  the  hybrid  Agl-Bi  aerosol  is  BbO;..  .V 
will  be  inferred  from  the  x-ray  analysis  it  is,  in  fact,  very  unlikely  the  incorporated 
form  of  bismuth  is  Bi203.  Hence,  the  resultant  aerosol  will  be  referred  to  as  Agl 
BiX  since  the  exact  form  is  unknown. 

As  mentioned  in  the  previous  chapter,  each  of  the  three  samples  were  identified 
from  an  x-ray  powder  analysis  as  iodargyrite,  the  hexagon  crystal  form  of  Agl. 
Neither  bismuth  nor  a  form  of  bismuth  was  detectable.  This  is  unusual  because, 
even  in  the  case  of  the  aerosol  generated  from  a  10  mole  °c  Bi I3  to  Agl.  the 
dominant  crystal  form  remains  hexagonal.  It  is  assumed  that  all  mole  ratios  of 
Bil3  result  in  iodargyrite.  Other  mixed  Agl  containing  aerosols  such  as  Agl-AgRr 
and  Agl-AgCl  have  a  face  centered  cubic  structure  (Vonnegut  and  Chessin,  1971; 
and  Palanisamy,  et  al..  1986).  Since  the  form  of  bismuth  does  not  change  the  type 
of  crystal  structure,  it  must  not  exist  as  BijtV  Bi203  has  its  own  crystal  structure, 
namely  monoclinic  with  a  pseudoorthorhombic  structure  type.  The  existence  of 
B12O3  would  preclude  the  formation  of  the  close  packed  hexagon  structure  of  Agl. 
Instead  the  form  of  bismuth  is  fully  incorporated  into  the  Agl  crystal,  and  this 
form  is  probably  inseparable  from  Agl  as  an  independent  molecule  or  crystal,  such 
as  Bi20„  or  BijO,,!,  (Personal  communication  with  Dr.  Oren  Anderson). 

The  changes  bismuth  imposes  on  the  crystal  structure  of  Agl  will  be  seen 
later.  First,  it  is  important  to  relate  the  d-values  of  Table  4.1  to  the  iodargyrite 
crystal  structure  (figure  5.1).  From  simple  geometry,  the  lattice  constants  along 
the  “c”  and  “a"  axis  are  determined  and  given  in  Table  5.1.  The  published  lattice 
constants  of  Agl  and  ice  are  given  for  comparison.  It  is  immediately  obvious 
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Table  5.1:  Lattice  Constants  from  X-ray  Analysis 


Sample 

“a"  Lattice  Constant 

“c"  Lattice  Constant 

Agl  aerosol 

4.616  ±  .002 

7.447  ±  .006 

AgI-0.5  mole  (x  Bi I3 

4.583  ±  .002 

7.493  ±  .006 

AgI-10.0  mole  9c  Bi I3 

4.622  *  .003 

7.567  r  .006 

*AgI 

7.5005  -  .003 

*Ice 

■ 

1 

7.367  -  .001 

*From  CRC  Handbook  of  Chemistry  and  Physics  standard  deviation  of  errors  (<* ) 
=  0.0019 


the  closest  match  to  ice  is  the  optimum  ratio  Agl-BiX  aerosol.  The  difference 
in  lattice  constants  between  the  Agl  aerosol  and  the  CRC  published  values  for 
Agl  lie  in  the  method  of  Agl  preparation.  That  is,  Agl  is  somewhat  unstable  and 
the  method  of  preparation  or  collection  (i.e.  precipitate,  quick-quench,  natural 
or  pyrotechnic)  may  greatly  affect  the  structure.  For  example,  the  preparation 
methodology  may  contribute  as  much  as  one  percent  oxygen  or  carbon  into  the 
final  Agl  structure  (personal  communication  with  Dr.  Dale  Winder).  But,  the 
most  important  conclusion  from  table  5.1  and  the  x-ray  analysis  is  that  the 
optimum  ratio  Agl-BiX  aerosol  has  the  closest  hexagonal  match  to  ice  ever 
reported.  Other  aerosols  have  shown  a  closer  match  to  ice  along  a  single  plane  such 
as  the  AgI-50  mole  9t  AgBr  aerosol  (Chessin  and  Yonnegut,  1971)  or  Metaldehyde 
(the  110  plane)  (Fukuta,  19G3).  However,  none  had  the  hexagonal  symmetry  more 
similar  to  ice. 

Another  important  point  to  draw  from  Table  5.1  are  the  differences  between 
lattice  constants  of  the  three  samples  with  ice.  These  differences  are  statistical!) 
significant  as  can  be  seen  from  each  constant's  error.  Table  5.2  emphasizes  these 
structural  differences.  The  disregistry  with  ice  is  the  percent  difference  between 
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Table  5.2:  Correlation  of  Effectiveness  with  Percent  Disregistry  with  Ice 


Sample 

fc  Disregistry 

Effectiveness 

a 

c 

(ice  crystals  g  nude 

Agl 

2.10 

2.58 

5.6  x  10'3 

Agl:  0.5  mole  /t 

Bi  I3 

1.34 

1.72 

4.4  x  10M 

Agl:  10  mole  °c 

Bi  I3 

2.23 

2.72 

4 ,0x  10,: 

Ice 

0 

0 

- — 

the  lattice  constants  of  each  sample  with  ice.  Again,  the  optimum  ratio  Agl- 
BiX  aerosol  has  the  better  match  with  ice.  One  can  only  speculate  as  to  how 
BiX  brings  the  Agl-BiX  crystal  into  a  closer  match  with  ice.  In  the  possible 
substitution  of  oxygen  and  bismuth  for  iodine  and  silver,  there  is  a  decrease  in 
the  ionic  radii  for  each  respective  substitution.  The  incorporation  of  a  foreign 
ion(s)  would  necessarily  cause  a  strain  in  the  lattice.  This  “lattice  distortion” 
will  be  smaller  or  larger  than  the  original  lattice  depending  on  the  ionic  radii  of 
the  foreign  ion(s)  (Palanisamy,  et  ah.  1986).  Agl-AgBr  and  Agl-AgCl  both  show 
such  a  “lattice  distortion”  in  substituting  smaller  ionic  radii  bromine  and  chlorine 
for  iodine  (Vonnegut  and  Chessin,  1971;  Palanisamy,  et  ah.  1986).  Figure  5.2 
gives  a  plausible  Agl-BiX  crystal.  The  left  portion  of  Figure  5.2  shows  a  model  of 
Agl  while  the  right  portion  shows  the  strained  Agl  with  a  bismuth  substituted  for 
one  silver  and  two  oxygens  substituted  for  two  iodines.  0“  and  BP3  have  ionic 
radii  of  1.40  A  and  1.20  A.  While  I-  and  Ag4-  have  ionic  radii  of  2.16  A  and  1.26 
A,  respectively.  The  result  is  a  smaller  crystal.  The  relative  size  of  each  crystal 
may  be  compared  ideally  by  overlaying  the  bottom  structure  on  top  of  the  upper 
structure  of  figure  5.2.  It  is  hypothesized  that  this  modified  Agl  crystal,  by  being 
a  closer  match  to  ice,  imposes  less  strain  during  the  epitaxial  growth  of  ice. 

In  relating  this  model  with  Table  5.2  and  Figure  4.1 .  as  the  amount  of  bismuth 
is  decreased  from  0.5  mole  percent,  a  fewer  number  of  Agl  aerosol  particles  art 
shrunk.  The  effectiveness,  in  turn,  approaches  the  effectiveness  value  of  Agl.  As 
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the  amount  of  bismuth  is  increased  from  the  optimum  A  g  I  -  B  i  1 3  ratio,  too  much 
bismuth  disrupts  or  strains  the  Agl  crystal  moving  the  structure  away  from  a 
match  with  ice.  Also,  excess  bismuth,  as  BijOs,  swamps  the  aerosol  population  a1' 
a  competing  aerosol  particle  with  poor  ice  nucleating  abilities. 

While  a  direct  mathematical  correlation  between  disregistry  and  effectiveness 
seems  obvious,  three  data  points  would  hardly  make  such  a  correlation  significant. 
This  is  an  area  where  further  research  could  be  done.  Other  evidence  and  specu¬ 
lation  as  to  the  nature  of  the  Agl-BiX  aerosol  will  be  discussed  with  reference  to 
heterogeneous  nucleation  theory. 

5.4  Heterogeneous  Nucleation  Theory 

Much  of  the  credibility  for  the  “closest  fit”  argument  is  the  theory  of  Turnbull 
and  Vonnegut  (1952)  which  relates  an  aerosol’s  disregistry  with  ice  to  its  threshold 
temperature  of  ice  nucleation.  In  an  attempt  to  relate  the  data  of  this  study  with 
Turnbull  and  Vonnegut’s  epitaxial  disregistry  theory,  a  bridge  between  the  theory 
and  experiment  was  sought.  Unfortunately,  effectiveness  data  is  too  inaccurate 
to  be  able  to  relate  to  the  minimum  temperature  of  ice  nucleation.  However,  the 
Agl-BiX  aerosol  clearly  owes  much  of  its  ice  nucleation  ability  to  its  close  match 
with  ice.  A  complete  x-ray  analysis  of  many  different  ratios  of  Agl  to  Bil3  could 
establish  this  relationship  in  much  greater  detail.  The  only  remaining  question 
may  relate  to  what  chemical  or  other  changes  might  have  occurred  during  the 
‘lattice  disruption’.  This  is  examined  in  the  next  paragraph. 

Another  aspect  of  heterogeneous  nucleation  theory  is  the  active  site  theory. 
The  data  of  this  study  presents  a  less  quantitative  picture  for  the  study  of  the 
active  site  theory.  However,  the  small  amount  of  BiX  incorporated  into  the  Agl 
lattice  mechanism.  A  small  ratio  of  a  compound  could  improve  the  nucleating 
ability  of  Agl  by  altering  or  augmenting  sites  to  make  them  more  hygroscopic 
(Corrin,  et  ah.  1967).  If  the  form  of  bismuth  is  an  oxide,  this  might  better 
facilitate  water  cluster  formation  (Zettlemoyer,  et  ah.  1963).  This  theory  also 
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better  explain  more  fully  the  quick  drop  off  of  effectiveness  with  ratios  of  Bil3 
greater  than  0.5  mole  percent.  Too  many  hygroscopic  sites  (i.e.  oxide  surfaces) 
might  create  a  water  monolayer  which  would  prevent  the  formation  of  water 
clusters.  Water  clusters  are  thought  to  be  a  possible  precursor  to  freezing.  Too 
many  hygroscopic  sites  might  also  favor  a  slower  deposition  mechanism  versus  a 
faster  contact  freezing  mechanism  (Edwards  and  Evans,  196S).  Figure  5.2  which 
served  as  a  plausible  structure  for  the  epitaxial  theory  also  serves  as  a  plausible 
model  for  the  active  site  scenario.  The  model  of  figure  5.2  exposes  two  oxides 
on  perpendicular  faces  of  the  iodargyrite  crystal.  But,  just  as  one  model  will 


explain  both  epitaxial  and  active  site  theories,  it  is  impossible  to  argue  totally  in 
favor  of  one  theory  at  the  expense  of  the  other.  As  the  crystal  structure  of  Agl 
shrinks  to  become  a  closer  match  with  ice,  sites  are  modified  and  brought  into  a 
better  juxtaposition  for  active  site  nucleation.  Thus,  both  theories  of  nucleation 
catalysis  are  probably  operating  in  concert. 


Chapter  6:  CONCLUSION 


6.1  Future  Research 

This  study  was  very  preliminary  in  its  examination  of  Agl  and  Bil3.  Other 
areas  were  beyond  its  intended  scope.  These  research  areas  and  several  questions 
remain  which  may  lead  to  future  research  efforts.  These  are  as  follows: 

1.  Particle  size  determination  was  not  done  in  this  study,  but  it  is  understood 
that  minor  differences  in  the  flame  and  generation  rate  compared  to  earlier 
studies  might  alter  the  distribution  of  particle  size.  This  could  in  turn, 
alter  the  relative  proportion  of  deposition  and  freezing  mechanisms  of  ice 
nucleation.  In  addition,  it  is  not  known  whether  Bil3  had  any  effect  on  the 
aerosol  size  distribution. 

2.  A  full  x-ray  analysis  is  needed  in  order  to  make  a  quantitative  study  of 
effectiveness  versus  lattice  disregistry. 

3.  A  complete  chemical  analysis,  working  along  with  the  x-ray  analysis  to 
completely  identify  and  characterize  the  Agl-BiX  aerosol  is  needed. 

4.  The  development  of  kinetic  theory  to  relate  threshold  temperature  to  an 
aerosol’s  rate  and  mechanism.  This  might  greatly  simply  future  studies  in 
characterizing  an  aerosol. 

5.  Improvement  of  the  laboratory  model  (such  as  correcting  boundary  effects) 
is  needed  before  a  proper  kinetic  theory  can  be  fine  tuned.  All  laboratory 
experiments  should  be  modeled  to  more  closely  approach  real  atmospheric 
conditions. 
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6.2  Summary 

One  of  the  original  aims  in  this  study  was  to  show  the  synergistic  relationship 
between  two  ice  nucleating  aerosols.  While  the  effectiveness  of  the  resultant 
aerosol,  Agl-BiX,  is  greater  than  the  effectiveness  of  Agl  or  Bi203,  the  synergistic 
relationship  obtained  was  not  quite  as  expected.  The  optimized  ratio  of  Bil3 
to  Agl  in  solution  was  determined  to  be  0.5  mole  with  an  effectiveness  an 
order  of  magnitude  greater  than  Agl  alone  at  -10CC.  Bil3,  which  had  previously 
been  identified  as  an  effective  nucleant,  oxidized  upon  generation  to  form  a  Bi203 
aerosol  which  was  ineffective  as  an  ice  nucleant.  The  rate  of  nucleation  with  Agl- 
BiX  did  not  differ  significantly  from  Agl,  and  contact  nucleation  was  identified 
as  the  ice  nucleation  mechanism.  However,  the  dramatic  rise  in  effectiveness  as 
the  ratio  of  Bil3  to  Agl  approached  0.5  mole  %  from  either  direction  (figure  4.1) 
prompted  an  investigation  into  the  reason  for  this  behavior. 

Previous  work  with  Agl  has  shown  the  unstable  nature  of  iodargyrite  (the 
hexagonal  crystal  form  of  Agl)  in  favor  of  the  face  centered  cubic  form  of  Agl 
when  combined  to  produce  Agl  containing  complex  aerosols  (such  as  Agl-AgBr 
and  Agl-AgCl).  However,  the  x-rav  powder  pattern  indicated  iodargyrite  even 
when  the  Bil3  ratio  was  as  large  as  10  mole  %.  It  was  assumed  that  the  hexagon 
form  of  Agl  was  maintained  at  all  mole  ratios  of  Bil3  tested  for  effectiveness. 

In  relating  the  laboratory  characterization  of  Agl-BiX  with  the  x-ray  powder 
analysis,  the  correlation  seemed  obvious.  The  optimum  ratio  Agl-BiX  aerosol 
had  the  closest  match  with  ice  along  both  the  “a”  and  “c"  axis.  The  “c"  and  ‘’a” 
axis  ice  disregistry  distinguished  the  Agl-BiX  aerosol  produced  via  combustion  of 
a  2/c  AgI-0.5  mole  %  BiI3-NH4I-acetone-water  solution  as  the  closest  hexagonal 
match  to  ice  ever  reported. 

Bismuth's  role  in  the  iodargyrite  crystal  is  hypothesized  to  create  a  ulatiic< 
disruption” .  In  shrinking  the  Agl  crystal,  the  larger  ionic  radii  of  silver  and  iodide 


are  substituted  for  smaller  ionic  radii  of  bismuth  and  oxygen.  Thus  less  strain  is 
encountered  in  the  epitaxial  growth  of  ice. 

The  correlation  of  effectiveness  and  disregistry  supports  the  closest  fit  or 
epitaxial  theory  of  heterogeneous  nucleation.  Clearly,  the  form  of  Bismuth  is 
repsonsible  for  the  change  in  the  Agl  lattice  to  more  closely  match  the  ice  crystal 
structure.  In  addition,  the  small  ratio  of  BiX  in  Agl-BiX  and  the  likelihood 
that  “X*  is  an  oxide  or  oxyiodide  supports  the  active  site  contention  of  nucle¬ 
ation.  Thus,  it  appears  heterogeneous  nucleation  functions  as  a  concert  of  subtle 
structural  and  surface  chemical  responses. 

Finally  the  implications  to  the  field  of  weather  modification  are  enormous. 
The  Agl-BiX  aerosol  can  be  used  in  all  ways  as  Agl  is  presently  being  used. 
More  important,  this  study  shows  how  a  very  small  and  specific  amount  of  Bil3 
can  significantly  improve  the  ice  nucleating  effectiveness  of  combustion  generated 
Agl.  Possibly  cheaper  contaminants  than  Bil3  could  be  chosen,  making  combustion 
generation  Agl  cloud  seeding  programs  even  more  cost  effective. 
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